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Using second-order gauge-invariant perturbation theory, a self-consistent framework describing 
the non-linear coupling between gravitational waves and a large-scale homogeneous magnetic field 
is presented. It is shown how this coupling may be used to amplify seed magnetic fields to strengths 
needed to support the galactic dynamo. In situations where the gravitational wave background is de¬ 
scribed by an ‘almost’ Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmology we find that the 
magnitude of the original magnetic field is amplified by an amount proportional to the magnitude 
of the gravitational wave induced shear anisotropy and the square of the field’s initial co-moving 
scale. We apply this mechanism to the case where the seed held and gravitational wave background 
are produced during inflation and find that the magnitude of the gravitational boost depends sig¬ 
nificantly on the manner in which the estimate of the shear anisotropy at the end of inflation is 
calculated. Assuming a seed held of 10 -34 G spanning a comoving scale of about 10 kpc today, 
the shear anisotropy at the end of inflation must be at least as large as 10 _4 ° in order to obtain a 
generated magnetic held of the same order of magnitude as the original seed. Moreover, contrasting 
the weak held approximation to our gauge-invariant approach, we hnd that while both methods 
agree in the limit of high conductivity, their corresponding solutions are otherwise only compatible 
in the limit of infinitely long-wavelength gravitational waves. 

PACS numbers: 98.80.Hw, 04.30.-w, 98.80.Cq 


I. INTRODUCTION 

The origin of cosmological magnetic fields that are prevalent throughout galaxies clusters, disk and spiral galaxies 
and high-redshift condensations has generated much debate in recent years, with the majority of this work being 
focused on providing mechanisms that generate these galactic fields on large scales (see [ij, 2] and references therein). 
The candidate mechanisms are diverse, often depending on the required seed field strengths. It has been suggested 
that the fields observed today could be a result of the amplification of a relatively large seed field through protogalactic 
collapse at the onset of structure formation Q. As the gas collapses to current measured densities, the flux lines of 
the frozen-in cosmological magnetic field get compressed, inducing adiabatic amplification. Another popular mech¬ 
anism, which requires a relatively weaker pre-existing seed field, is amplification via the galactic dynamo by means 
of parametric resonance Q. The combined effect of differential rotation across the disk and the cyclonic turbulent 
motions of the ionized gas is believed to lead to the exponential amplification of a smaller primordial held until the 
back-reaction of the plasma opposes further growth. Although the dynamo mechanism is strongly supported by the 
close correlation between the observed structure of the galactic fields and the spiral pattern of galaxies, there is some 
argument over its efficiency and hence the amount of amplification that can occur through this process. The major 
problem with all of these mechanisms is that they assume the presence of a pre-existing seed field whose origin is 
still to be established. A further idea relies on turbulence (disrupted flow) and shocks, which occur during the stages 
of structure formation, inducing weaker magnetic fields via battery-type mechanisms, which operate as a result of 
large-scale misalignments of gradients in electron number density and pressure (or temperature) ||. 

There have been numerous attempts to generate early, pre-recombination, magnetic fields with strengths suitable 
to support and maintain the dynamo by exploiting the different out-of-equilibrium epochs that are believed to have 
taken place between the end of the inflationary era and decoupling jjj]. These fields are facilitated by currents that 
arise from local charge separation generated by vortical velocity fields prevalent in the early plasma (cf. also 0 )- 

One problem with the above mechanisms is that they are casual in nature so the scales over which the fields are 
coherent cannot exceed the particle horizon during that epoch. Given that such phase transitions took place at very 
early times, where the comoving horizon size was small, tight constraints must be placed on the coherence length of 
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these magnetic fields. However, pre big bang models based on string theory Q, in which vacuum fluctuations of the 
magnetic field are amplified by the dilaton field, predict super-horizon fields. 

Inflation has long been suggested as a solution to the causality problem, since it naturally achieves correlations on 
superhorizon scales, however adjustments to the standard inflationary models need to be made since magnetic fields 
surviving this epoch are small on account of the inability of vector fields to couple gravitationally to the conformally 
flat metric resulting from the exponentially fast expansion. A way around this obstacle is by breaking the conformal 
invariance of electromagnetism since this alters the way the underlying gauge fields couple to gravity. There are many 
ways of doing this which explains the variety of the proposed mechanisms in the literature Q. Such inflationary 
scenarios have not been without critique, though Ha¬ 
lt has also been proposed that inflation is followed by a period of preheating in which the parametric resonance of 
the causal oscillations of the inflaton field and the accompanying perturbations can lead to amplification on super¬ 
horizon scales 0 - Other authors have advocated the breakdown of Lorentz invariance either in the context of st ring 
theory and non-commutative varying speed of light theories, or due to the dynamics of large extra dimensions [121 . 
The success of these proposals, however, is usually achieved at the expense of simplicity. 

In order for these proposed mechanisms to be viable, they must, in addition, produce seed fields that satisfy the 
criteria for the subsequent amplification processes to work. To be a candidate seed field for the galactic dynamo, the 
induced field must exceed a minimum coherence scale in order to prevent the destabilization of the dynamo action. 
The time scale over which the amplification takes place also dictates a minimum field strength, for example in the case 
of a dark-energy dominated Universe we obtain B ~ 10 -34 G on a coherence scale of 10 kpc. Davis et al. 0 proposed 
an inflationary mechanism that exploits the natural coupling between the Z-boson and the gravitational background. 
Unfortunately, the fields produced only just fall within dynamo limits in the case of a dark energy dominated Universe. 
Recently, the production of a magnetic seed field due to the rotational velocity of ions and electrons, caused by the 
nonlinear evolution of primordial density perturbations in the cosmic plasma during pre-recombination radiation and 
matter eras, was investigated in 0 and a rms amplitude B ss 10 23 (A/Mpc) 2 G at recombination on comoving 
scales A > 1 Mpc was reported. 

In this paper, we offer an alternative mechanism that looks at the interaction of a pre-existing field, such as 
the one proposed by Davis et al., with a gravitational wave (GW) spectrum which accompanies most inflationary 
scenarios. This builds on earlier work by Tsagas et al. 0 in which this idea was first introduced within the weak 
field approximation. Our aim is to investigate whether this interaction can produce a sufficiently large amplification 
of a seed field present at the end of inflation to meet the above mentioned requirements for the dynamo to work. 

The issue of how to deal with the coupling between gravitational waves and the seed magnetic field is rather subtle. 
A commonly used approximation in the literature is to assume that the magnetic field is weak and that its contribution 
to the energy-momentum tensor is such that it does not disturb the isotropy of the FLRW background 16]. This is 
done by assuming that the energy density of the magnetic field B a is much less than the matter energy density: B 2 -C p 
and that its anisotropic pressure is negligible: ir a b = — B <a Bb> « 0 12Pl | . The problem with this approximation is that 
it is not gauge-invariant in a strict mathematical sense, so one can therefore not guarantee that, when calculating the 
magnetic field which arises through its coupling with linear perturbations of FLRW (such as gravitational waves), it 
leads to physically meaningful results. In order to solve this problem we develop a self-consistent framework based on 
second-order perturbation theory, employing the methods initiated by recent work of Clarkson 0 and Clarkson et 
al. 0 - Here the seed magnetic field is treated as a on average homogeneous linear perturbation of the background 
FLRW model and couplings to gravitational degrees of freedom that arise when perturbing the background are taken 
to be second order in the perturbation theory. Adopting this approach allows us to write Maxwell’s equations in 
a way that makes them manifestly gauge-invariant to second order with interaction terms that clearly describe the 
modes induced by the gravity wave-magnetic field interaction. The restriction to a homogeneous seed field leads to 
simplification on the technical level but still encapsulates the main features of the gravito-magnetic interaction. The 
implementation of an inhomogeneous seed is reserved to a future article. 

The results show that, in the presence of gravitational radiation, the magnitude of the magnetic field is amplified 
proportionally to the shear distortion caused by the propagating waves. Once the amplification is saturated, the 
magnetic field dissipates adiabatically as usual. The gravitational boost is also proportional to the square of the 
field’s original scale, which suggests that the proposed mechanism could lead to significant amplification in the case 
of large scale magnetic fields. Indeed, when applied to fields of roughly lO -34 G spanning a comoving scale of about 
10 kpc today (see for example the fields produced in [u|) , the mechanism leads to an amplification of up to 13 orders 
of magnitude (depending on the calculation of the shear distortion), bringing these magnetic fields well within the 
galactic dynamo requirements, without the need for extra amplification during reheating. We thus qualitatively and 
quantitatively rediscover in a gauge-invariant fashion the main results reported in 0 . 

In order to contrast the two different approaches in detail, we compare our solutions with the corresponding 
solutions obtained using the weak field approximation lid and find that while both methods agree in the limit of high 
conductivity, their corresponding solutions are otherwise only compatible in the limit of infinitely long-wavelength 
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gravitational waves when merely the dominant contribution is considered. 

The units employed in this paper are c = h = 1 and n = 87 tG = 1, the exception being section El where natural 
units are used. 


II. PERTURBATION SCHEME 

If we wish to study the interaction between gravitational waves and a magnetic field in a cosmological setting, we 
immediately face a second-order problem in perturbation theory because both the magnetic field as well as GW are 
absent in the exact FLRW background, and may thus be individually regarded as first order perturbations. Using 
the 1+3 covariant approach |lfl| . we therefore develop a two parameter expansion in two smallness parameters: cb 
represents the magnitude of a homogeneous magnetic field and e g represents the magnitude of the GW. The magnitude 
of the interaction GW x magnetic field is of order 0{tBtg) as is the magnitude of the in such a manner generated 
electromagnetic fields. However, at second-order level, only terms of order 0{tB^g) are kept while terms of order 
0(e 2 ) and 0(e 2 B ) are discarded. In fact, when dealing with the gravito-magnetic interaction, these discarded terms 
would always appear multiplied by a first-order quantity and are thus irrelevant for our considerations. 

Whence, the perturbation spacetimes are divided up and denoted in the following way: 

• B = Exact FLRW as background spacetime, C>(e 0 ); 

• T\ = Exact FLRW perturbed by a homogeneous magnetic field whose energy density and curvature are neglected, 

O(eu); 

• Ti = Exact FLRW with gravitational perturbations 0(e g ); 

• S = T\ + Ti allows for inclusion of interactions terms of order Oiestg)- 

We will generally refer to terms of order 0(cb) and 0(e g ) appearing in T as ‘first-order’ and to variables of mixed 
order 0{es^g) appearing in S as ‘second-order’. 

It should be noticed that the absence of an electric field in T\ does not necessarily imply that there is no electric 
field at all but rather that the electric field is perturbatively smaller than the magnetic field. This is in accordance with 
the standard assumption that the very early Universe was a good conductor (see, for example, pp| for an example 
of how this works). The inclusion of an electric field in T\ is possible, in principle, but would require to alter the 
perturbation scheme because then interactions between gravitational waves and the electric field needed to be taken 
into account as well. However, a more realistic way of describing the interaction between gravitational waves and 
electromagnetic fields should employ a multifluid description [?j], which allows for modelling the currents, but that is 
beyond the scope of the present paper. 

Having outlaid the different stages we turn to review the concomitant equations. We keep them as general as 
possible, which will allow us to illuminate the effects of spatial geometry, cosmological constant A and equation of 
state for the matter on the interaction. We limit ourselves to the irrotational case, that is, we require the vorticity 
Lo a b to vanish throughout. 


A. FLRW background 

The FLRW models are characterized by a perfect fluid matter tensor and the condition of everywhere-isotropy. 
Thus, relative to the congruence of fundamental observers with 4-velocity u a (u a u a = — 1), the kinematical variables 
have to be locally isotropic, which implies the vanishing of the 4-acceleration ii a = u b \/bU a , shear a a b = D < a Ub> 
and vorticity oj a b = D[ a Ub] (0 = u a = <J a b = u a b)- Furthermore, the models have to be not only conformally flat, that 
is, the electric and magnetic components of the Weyl tensor vanish (0 = E a b = H a b), but also spatially homogeneous 
implying the vanishing of the spatial gradients of the energy density pt, the pressure p and the expansion 0 = D a u a 
(0 = D a /i = D a 0 = D a p). As usual, the spatial derivative D a = h b V& is obtained by projection of the spacetime 
covariant derivative V a onto the 3-space (with metric h a b = g a b + u a Ub ) orthogonal to the observer’s worldline. As a 
consequence, the key background equations are the energy conservation equation 

/i + 0(/x + p) = 0 , (1) 


the Raychaudhuri equation 


0 — —^0 2 — \ (m d - 3 p) + A , 


(2) 
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and the Friedmann equation 


a = 5 02 


3 K 


where the constant K indicates the geometry of the spatial sections. 


( 3 ) 


B. First-order perturbations 

1. The homogeneous magnetic field B a 

We assume the magnetic field B a to be spatially homogeneous at first order (D a B b = 0) and thus consider the 
gradient of B a as well as the magnetic anisotropy n o b = —B <a B b> as being of second order. We presuppose that 
such a field was produced by some primordial process, which left a relic field on average homogeneous over a typical 
coherence length. Since there are no electric fields or charges in the T\ perturbation spacetime, the magnetic induction 
equation takes the form 


Pa = B <a> + |©-B a — 0 . 


( 4 ) 


As a result, the magnetic field scales as 



where a denotes the scale factor, e.g., 0 = 3 a/a = 3 H, where H denotes the inverse Hubble length. 


( 5 ) 


2. Gravitational waves 

Gravitational waves are covariantly described via transverse parts of the electric (E ab ) and magnetic (H ab ) Weyl 
components, which are PSTF tensors l2l| . The pure tensor modes are transverse, obtained by switching off scalar 


and vector modes (0 = D a /x = D a @ = D a p = cu a = u a ), which results in the constraints |3(j 

0 = D a cr ab = D a E ab = D a H ab = H ab - curl a ab ■ (6) 

The propagation equations for these tensor modes are simply 

d<ab> + g0Cr a b — E ab , (7) 

E <ab >+QE ab = curl (curl a ab ) - ^ (jt + p) cr ab , (8) 

together with the background equations for 0 and fj,. Since every FOGI tensor satisfies the linearized identity 

curl (curl T ab ) = —D 2 T ab + |D <a D' T b>c + (/z + A — ^0 2 ) T ab , (9) 

we see that the gravitational waves are completely determined by a closed wave equation for the shear, namely 

dab — D 2 (J ab + |0(j a 6 + (^0 2 + + |A) (Tab = 0 . (10) 


C. The interaction 


Maxwell’s equations govern the interaction between GW and magnetic fields. If we require charge neutrality and 
neglect currents as well as the back-reaction of induced second-order magnetic fields with the shear, we obtain 


E<a> + §0-Ea — CUrlH a , 

B<a> + §0-B a = (JabB b - curl E, 

D a E a = 0, 


( 11 ) 

( 12 ) 

(13) 

(14) 


D a B. 
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Observe that the EM fields have to be divergence-free at all orders due to neglecting vorticity effects. Moreover, the 
system is not gauge-invariant because it contains a mixture of second-order (E a , curl E a: curl B a ) and first-order terms 
(<7 a &), while B a now comprises the full magnetic field (the first-order contribution plus the induced field). The situation 
we are interested in is the interaction between the shear a a b and the first-order magnetic field, neglecting the back- 
reaction with the induced magnetic field. How does one then disentangle the different magnetic field perturbations in 
a consistent way? 

In special relativity, the standard procedure would be to use a power series expansion of the magnetic field, 

B a = cbB* + e g e B B a 2 + 0(e 2 g , e 2 B ) , (15) 

where the first-order field £?“ satisfies the magnetic induction equation 0 . Although insertion of this expansion 
into the above system yields only second-order terms, the procedure does not work in general relativity since the 
commutation relations for the various differential operators (cf. the appendix) can not be consistently satisfied. To 
illustrate this important point clearly, we consider the commutation relation between the (proper) time derivative 
and the spatial gradient applied to the magnetic field. It is evident that the case where the commutator relation is 
introduced after the expansion of B a , 


{D b B“) ± =e g e B (D b B*) ± 


e g e B 


D b BZ - \QD b B9, 


(16) 


does not agree with the case where the linearized identity for ( D a B b ) is substituted before using the power series 
expansion m- 


( D b B a ) ± 


D b B a - \QD b B a + H bd e dac B c + a d c D c B a 




D b B% - \QD b B“ 


e B H b e dac Bl 


(17) 


Here, _L denotes projection onto the fundamental observer’s rest space. This inconsistency can only be resolved if all 
interaction terms are zero. It is via the commutation relations that Weyl curvature is brought in through the back 
door which couples to the magnetic field and thus affects the interaction. It is this feature that renders the power 
series procedure faulty. 

The difficulty arises because the magnetic field B a is not gauge-invariant in S as it does not vanish in T\. We 
therefore need to define a new second order gauge-invariant (SOGI) variable which satisfactorily describes the effects 
that we wish to investigate. However, a look at Maxwell’s equations above reveals that (3 a = B <a> + | QB a is the 
sought SOGI variable which has to be used at second order instead of the magnetic field B a . We chose to describe the 
interaction in terms of the variable I a = <J a bB b . Hence, Maxwell’s equations can be written in truly gauge-invariant 
terms at second-order, namely 

E <a > + %OE a = curl B a , (18) 

(3 a + curl E a = I a . (19) 


Observe that the standard constraints 0 = T) a B a = D a E a , which hold at all orders, imply 

D°/3 q = D a I a = a ab B a B b = 0 , 


( 20 ) 


where the last equality is only true as long as spatial gradients of B a are regarded as second-order. Clearly, if the 
idealized assumption of infinite conductivity is made so that all electric fields vanish, Maxwell’sequations reduce to 
(3 a = la- In this specific case, once the solution for I a is known, the (not gauge-invariant) generated magnetic field 
measured by the fundamental observer can be obtained via a standard integration of (3 a . However, it is important to 
stress that (3 a is the fundamental variable, whose deviation from zero quantifies the evolution of the magnetic field at 
second-order in a truly gauge-invariant manner. 


III. WAVE EQUATIONS FOR THE MAIN VARIABLES 

Having written the key Maxwell’s equations as a system of differential equations of purely SOGI variables, we now 
turn to the derivation of wave equations for the electric and magnetic fields. In doing this we make no assumptions 
about the spatial geometry or the equation of state and also keep the cosmological constant; this has the advantage of 
allowing us to draw some conclusions about how these parameters influence the interaction between GW and magnetic 
fields. In particular, it will turn out that neglecting the current in Maxwell’s equations and at the same time requiring 
a homogeneous magnetic field at first-order level leads to consistent equations in spatially flat models only. 
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A. Wave equation for the interaction variable 

Let us first derive the wave equation for the interaction variable I a = <J a bB b . Even though the shear a ab belongs to 
J~2 and the magnetic field B a to the commutator relations do not lead to ambiguities for I a since they manifest 
themselves only at third-order in this case. In order to derive an evolution equation for I a , we need the auxiliary 
quantity J a = E a bB b . Then, using equations 0 , 0 , 0 and < 0 , we arrive at the system 

I<a> + f ©4 = - Ja , ( 21 ) 

j<a> + §©4 = — D“7a + [5 (/X — p) + A — ^©“] hi , (22) 

where we employed that spatial gradients of the magnetic field are second-order and thus D 2 / a = D 2 (a a bB b ^J = 
(D 2 ( 7 a f,) B b . Eliminating the auxiliary variable J a , the general closed wave equation for I a is found to be 

I<a> ~ D 2 Ia + 30/<Q> + [qy © 2 — — \p + 5 A] I a = 0 . (23) 

In the case of infinite conductivity, the solution to equation PI) instantly yields the solution of 0 a , from which the 
induced magnetic field measured by the fundamental observer might be obtained by integration. 


B. Wave equation for the electric field 

To derive the wave equation for the induced electric field, we first differentiate equation CHI and equate the result 
with the second-order identity 

(curl B a ) ± = -© curl B a + curl (3 a - H ab B b (24) 

to obtain 

E <a > + §©£«*> + [f © 2 - | (/r + 3 p) + |A] E a = curl 0 a - H ab B b . (25) 

Secondly, using equation PI) to substitute for curl 0 a above and the expansion 

curl (curl E a ) = -D 2 E a - [§0 2 - f (ji + A)] E a , (26) 

we find a forced wave equation for the induced electric field, namely 

E <a> - D 2 E a + l&E <a> + [|0 2 + 1 ( M - 3p) + |A] E a =K a , (27) 

where the forcing term I\ a = curl I a — H ab B b = e c d[ a D rr fc j c B b has no divergence. It is possible to show that the forcing 
term K ai as well as curl I a and H ab B b , respectively, can be found from the wave equation 

K <a > - D 2 A' q + f ek <a> + [f 0 2 - i (/r + 9p) + |A] I< a = 0 . (28) 

For example, the wave equation for curl I a follows by taking the curl of equation (El and using the expansion El , 
while the case H a bB b is similar to the derivation of the wave equation for the interaction term I a . 

It will be useful for later purposes to consider the electric field’s rotation. By taking the curl of equation E3 , we 
immediately arrive at 

(curl E a )_i - D 2 (curl E a ) + |© (curl E a ) ± 

+ [|© 2 + | (ft - 9 p) + |A] curl E a = curl K a . (29) 

Because curl ^ H a bB b ^ = — D 2 I a + [— y^© 2 + | (p + A)] I a holds, we note the interesting result 

curl K a = [jb© 2 - i (n + A)] I a . (30) 

That is, for a cosmological model with flat spatial sections we have curliv 0 = 0 and, therefore, the electric field’s 
rotation is not induced by the interaction between magnetic fields and GWs at second-order - the generated electric 
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field is curl-free. As a consequence, the interaction between magnetic field and GW produces the same magnetic field 
in a case of a spatially flat Universe as in the limit of high conductivity. 

However, upon closer inspection of the forcing term K a in equation m one discovers that this term is actually 
identically zero because of the identity fzl 

0 = e abc V b (D d A c d )-2V b e cd ^(v c A b] d ) , (31) 

which holds for any vector V a and tensor A a b = A <a b > perpendicular to the congruence u a . Thus, equation m 
implies that our chosen perturbative scheme is only consistent if the cosmological model is spatially flat (cf. also 
footnote B below). In essence, we see that the requirement of having a spatially homogeneous and thus a curl-free 
magnetic field at first-order can only be achieved when the Universe is spatially flat. Furthermore, the interaction 
between GW and a magnetic held generates in this particular case no electric Helds (at least to second order in the 
perturbation scheme). 


C. The generated magnetic field 

We have already pointed out that for spatially fiat models the generated magnetic field follows directly from the 
interaction variable since in this case we have p a = J a . For closed or open models, however, a wave equation for 
P a is needed to determine the induced magnetic field. The sought after equation may be obtained by applying the 
constraint equation m to equation l|29|) and substituting for curl I\ a via equation m- This leads to 

P<a> — D 2 Pa + f©/3<a> + [|© 2 + g (t 1 — $P) + §A] /3 a = 

I<a> — D 2 I a + + [y§0 2 + ~ §V + 7 f A] I a • (32) 

Observe that for models with flat spatial sections the lhs and rhs of the above equation become identical - in agreement 
with the comment following equation A slight simplification is achieved by employing equation (ESI yielding 

finally a forced wave equation for p a : 

P<a> ~ D 2 Pa + ^Qp<a> + [|© 2 + g (M — $P) + f A] Pa = 

- Ie/< a > - [§0 2 - § a* + 2 p- A)] I a . (33) 

It is evident that the variable I a and hence the gravitational waves source fluctuations in the magnetic field variable 
P a . Another way to derive equation consists of differentiating Maxwell’s equation ni twice, using equation ® 
to get rid off the curl U a -term and applying the corresponding commutation relations. This clearly demonstrates the 
consistency of our approximation scheme. 


IV. SOLUTIONS FOR FLAT UNIVERSES 

After having derived the fundamental equations governing the interaction between GWs and magnetic fields as well 
as the generated electromagnetic fields, we turn to the task of solving them. For the sake of simplicity, we investigate 
the solutions only for spatially flat models with zero cosmological constant A. We assume the matter to obey a 
barotropic equation of state, p = wp, with constant barotropic index w. 


A. A useful time variable 

The background equations dll subject to the assumptions stated above imply the following evolution equation 
for the scale factor: 


- + i(l + 3u0(-) =0. (34) 

a 2 \a J 

By integrating once and choosing initial conditions such that @o = B(to) = 3Ho f° r some arbitrary initial time to 
with H = a/a defining the Hubble radius, we obtain the following solution for the expansion 


ie = “ =_ l _. 

3 a 3(1 + w)(t — to) + 2 /Hq 


(35) 



Integrating once more, we find for the scale factor the solution 


a(t) = a 0 [§ H 0 (1 + w)(t- t 0 ) + l] 3 ( 1 +* u ) . (36) 

The introduction of a dimensionless time variable, r, defined for w ^ — 1 as 

T=§J?o(l + «0(i-io) + l , (37) 

will turn out to be extremely useful as it simplifies the integration of almost all the equations considered later 
irrespective of the barotropic index while taking the initial conditions explicitly into account as well. For example, 
the scale factor evolves simply as a = ao r 2 ^ 3 * 1+u ’^ and the Hubble radius is given by H = Hq/t. Moreover, r = 1 
corresponds to the initial time to- Note however that the r variable cannot be used in the de Sitter limit w —> — 1, 
where the scale factor becomes a(t) = aoexp (H 0 (t — to)). 


B. Generated magnetic field 


Since we are only considering Universes with flat spatial geometry, the induced magnetic field can be found by 
integrating over j3 a . To this end, it suffices to solve for the interaction variable I a . A standard harmonic decomposition 
is used to take care of the Laplacian operator. We expand the shear a a b = Ylk a ^Qab 3n P ure t ensor harmonics, 
where as usual Q<l b> = 0 and D 2 Q^ b = — (k 2 /a 2 )Q^ b hold. Moreover, each gravitational wave mode is associated 
with the physical wave length Aqw = 2iTa/k. Since the magnetic field in T\ obeys curlH a = 0, it follows that 
D 2 B a = —curl (curl B a ) = 0 and therefore that the expansion of the magnetic field B a = J2 n B^Qa in pure vector 

(solenoidal) harmonics reduces to B a = B^Qa , where B = B°(ao/a) 2 . This just means that the magnetic field B a 
is spatially constant, e.g., in agreement with the assumption of homogeneity |inj. Of course, the solenoidal harmonics 
also obey the relations = 0 and D 2 Q^ = —(n 2 /a 2 )Q^\ Perturbations in S are conveniently decomposed with 

the vector harmonics [32] Va ' ] = Q\ n ) , which are readily verified to fulfill the standard requirements = 0 

and D 2 Va £> = —(£ 2 /a 2 )Va f \ where the wavenumber £ satisfies £ 2 = ( k a +n a )(k a + n a ). Because the magnetic field in 
T\ has got only the zero mode in our investigation, the wavenumber t coincides with the wavenumber k of the shear. 

Using the unified time variable r and the harmonics explained above, we transform the wave equation m for the 
interaction variable I a into an ordinary differential equation: 


^(1 + w) 2 r; 


w 


27(1 +w) 

2 T 


+ 


aoHo 


T 3(l + m) 4- 


25 - Ibw 


2 T 2 


1(e) = 0 , 


(38) 


where a prime means differentiation with respect to r. Initial conditions are chosen as follows: 

1(e) (to) = &(k)(to)Bo , 


I(e)i T = 1 ) = B 0 


^ Cfc) (i) - 


3(1 + w) 




here, Bq is the initial amplitude of the first-order magnetic field and 

<j(k)(to) = 3/2Ho (1 + w) cr( fc) (l) 

was used. For every mode k we have initially <7(to) = oq and cr’{1) = a r 0 . 


(39) 

(40) 


(41) 


1. Infinite-wavelength limit 

In the infinite-wavelength limit (£ —> 0), the solution of equation (13811 is easily found to be 

//-jN _ 10 —5^310 

I (0 > ( T ) = Ci t + C 2 , 


(42) 
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where C\ and C 2 are constants of integration. If the initial conditions CUB PB are chosen, the corresponding 
integration constants are 


(-5 + 3 W )J ffl (l)-3(l+ W )/( £) (l) 
5 + 3w 

10I W (1) + 3(1 + w)I' w (l) 

5 + 3 w 


(43) 

(44) 


We remark that this solution is in agreement with the result obtained by multiplying the first-order magnetic field 
© with the infinite-wavelength solution of the shear equation m- Whence, the total magnetic field in the presence 
of infinite-wavelength GWs is 


B (0) (r) = B 0 t~^ rW 



Ci 2 

B 0 H 0 3(1 — w) 


( -1+™ 
T 3(1+w) 



C 2 1 
B 0 H 0 1 + 3w 


/ 2+6™ 

\ 

( r 3(l+u,) 



(45) 


where Bo is the magnitude of the first-order magnetic field interacting with the GW at initial time to and it is required 
for physical reasons that the induced magnetic field vanishes initially. We stress that the interaction always leads to 
an amplification of the magnetic field for any physically acceptable choice of equation of state because of the growing 
contribution in the second line of equation (TTIfll . 

Let us look at some important special cases. For the sake of simplicity, we take 1^(1) = 0 for granted. In the 

matter-dominated era, where the matter is accurately described as dust, w = 0 and a = ao t 2 / 3 , this yields for the 
magnetic field mode 


<„<«)= a, (%y 

whereas for a radiation-dominated era, where 



1 + 

w = 

1 + 


2 CTo_ f / Oo \ 3 / 2 \ 2a 0 f a V 

3H 0 \\a) J + H 0 \a 0 iJ ’ 

1/3 and a = oqt 1 / 2 , the magnetic field mode is 



(46) 


(47) 


It follows that in the infinite-wavelength limit the amplification depends mainly on the scale factor and the magnitude 
of the initial GW distortion relative to the Hubble parameter (<j/H) 0 - 


2. General case with l ^ 0 


The general solution to the interaction equation ifTSll is: 


I{€)(t) = 


Di Ji 


3w + 5 l 

9 ( 1 -I- Hr, Hr, 




1 + 3™ 

_ r 3(1+™) 


- Do .To 


( 3^ + 5 £ 


l+3u 


(48) 


where D 1 , D 2 are integration constants and Ji, J 2 denote Bessel functions of the first and second kind, respectively. 
Observe that in the limit of infinite wavelengths, £ —> 0, the solution (PI) is recovered. The generated magnetic field 
relative to the observer moving with 4-velocity u a can be calculated from the solution (HQ) analytically for every 
barotropic parameter w. We will state here only the total magnetic held solution in the case of dust and radiation, 
respectively. For dust, where w = 0 and a = qqt 2 / 3 , the full magnetic held is 


B 


W 

Dust 


(a) = Bo (^)‘ 


1 + 


3 f Agw Y 

47T 2 \ Ah / ( 


Ho 


2//i 


+ 0(a~ 1 ) 


while for radiation, where w = 1/3 and a = aoT 1 / 2 , the total magnetic held modes obey 






(49) 


(50) 


Here, we introduced the gravitational wavelength Aqw = 27t a/k and the Hubble length Ah = 1 /H. The un-displayed 
remainders 0(a~ 1 ) in the expressions above contain oscillating functions which decay at least as fast as the inverse 
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scale factor a 1 . Note that when the infinite-wavelength limit of the full solutions above is taken, the findings 
and m are rediscovered. The results (TTT71) m clearly show how the generated magnetic field depends on the initial 
conditions and that the late time behaviour is almost identical for both dust and radiation. It should be noted that the 
interaction can only be effective if the wavelength of the GW matches the size of the magnetic field region, Agw ~ Ag: 
in the case of Agw Ag the magnetic field cannot be physically affected by the GW, while for Agw “C Ag the effect 
becomes negligible due to its quadratic dependence on Agw- If we divide the findings m ct through the energy 
density of the background radiation, the dominant contribution can be summarized as follows, 



where the wavenumber indices have been suppressed and a' 0 = 0 was assumed. At late times, a significant amplification 
of the original magnetic field can be achieved for super-horizon gravitational waves. Note that a result almost identical 
to 15111 was obtained in 0|, wherein the factor 1/10 is replaced by 10 instead. However, our result holds for any finite 
gravitational wavelength, Agw ~ Ag, while the result in 0 assumes Ah Agw- Moreover, 0 used somewhat 
contrived initial conditions leading to an abrupt amplification of the field whereas we chose initial conditions such 
that there is no generated magnetic field present when the interaction kicks in at the end of inflation. 


V. APPLICATION 


In order to estimate the amplification of the seed field due to the interaction with GWs we reproduce the analysis 
presented in m using the same parameter values. We find it convenient to adopt natural units in this section. 

Given that the evolution of the (spatially flat) Universe is dominated by a dark energy component such as a 
cosmological constant or quintessence, the minimum seed required for the dynamo mechanism to work is of the order 
of 10 _3 ° G at the time of completed galaxy formation and coherent on a scale at least as large as the largest turbulent 
eddy, roughly ~ 100 pc 0|. Such a collapsed magnetic field corresponds to a field B of ~ 10 -34 G with coherence 
length Ag ~ 10 kpc on a comoving scale if the field remains frozen into the cosmic plasma from the epoch of radiation 
decoupling to galaxy formation. Its field strength compared to the energy density of the background radiation, /i 7 , 

gives rise to the ratio B /[u]' 1 ~ 10” 29 , which stays constant as long as the magnetic flux is conserved and the magnetic 
field is frozen into the cosmic medium. 

During inflation, the Hubble parameter H remains constant and is taken to be H ~ 10 13 GeV 0 . The scale of 
the magnetic field therefore implies Ag/An ~ 10 20 at the end of inflation. A general prediction of all inflationary 
scenarios is the production of large scale gravitational waves whose energy density per wavelength is roughly fulll 


Ww^m^(-1) 2 (-V . (52) 

\Agw/ \mpij 

Here, Aqw denotes the wavelength of the GW and mp\ the Planck mass (see, for example, HU). The total energy 
density of the gravity waves expressed in terms of the shear is (see footnote 4 in |25| for a neat discussion) 


m pi 

^GW = TV- O 'ab cr 
107r 


ab 


which implies an induced shear anisotropy El 



(53) 


(54) 


where the zero suffix indicates the end of the inflationary epoch. Typical inflationary models predict H/mp\ ~ 1CU 6 , 
which lies comfortably within the bound H/mp\ < 1CD 5 stemming from the quadrupole anisotropy of the CMB. 

The interaction of such a primordial magnetic field with GWs produced by inflation leads to a substantial amplifi¬ 
cation of the former. Resorting to our result m and applying (151) . we find for the magnetic field El 


B 






(55) 


Substituting (Ag/An)o ~ 10 20 and H/mp\ ~ 10 6 into the above expression, we obtain that GWs amplify the original 
magnetic field as much as 13 orders of magnitude. This mechanism thus brings an inflationary seed such as in El 
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up to ~ 10 21 G, which is comfortably within the requirements of the galactic dynamo mechanism (23j]. In Universes 
with zero cosmological constant, the minimum seed for the dynamo has to be raised from ~ 10 -30 G to ~ 10 -23 G [26^ . 


However, the use of and innt to find the shear anisotropy is problematic in the sense that P i holds 

strictly speaking only on scales up to the horizon; the energy stored in superhorizon modes cannot be measured by 
local observers. Once the initially superhorizon GWs re-enter the observer’s horizon, they contribute to the measured 
energy density. Therefore, the correct procedure is to use the value of the shear anisotropy at horizon crossing, 
(ct/H) H q, and scale that value back to the end of inflation using its evolution equation. During the radiation 
dominated era, the gravitational wave length varies with the scale factor (i.e. Aqw ~ r 1 / 2 ), while the horizon scales 
as the inverse of the square of the scale factor (i.e. Ah = H~ x ~ r). This gives 


f Agw\ _ / Agw\ _i/2 

Uh/ Uh /o 


(56) 


where r = 1 corresponds to the end of inflation after reheating. At the point in time where the gravitational 
wave crosses back into the horizon, its physical wavelength equals the Hubble scale Ah, that is (Agw/Ah)hc = 1- 
Substituting for the values from above, we obtain the crossing time thc = 10 40 , which hitherto leads to the Hubble 
parameter Hhc ~ 10 _27 GeV at horizon cr ossing. Du ring the radiation era, the temperature is proportional to the 
square of the Hubble parameter 0, T ~ ^/mpL-ff/10, yielding a temperature of Thc ~ 10 -5 GeV at horizon crossing. 
Since 1 GeV ~ 10 13 K, the actual temperature is ~ 10 8 K, confirming that the superhorizon mode under consideration 
indeed crosses back into the horizon during the radiation-dominated epoch. Whence, combining Pi and m gives 
a shear anisotropy 


Shc = 



(57) 


at the time of horizon crossing. 

Assuming a flat model with no cosmological constant, the evolution of the shear anisotropy, £ = <j/H , can be 
obtained by solving for the shear modes a^ :) from equation 1101 ) and noting that H ~ r _1 . The result is 


-1 + 9uj 

S (fc)( T ) = T®' 14 ”’ 


AJi 


3 w + 5 


2 (1 + 3 w) ’ ao-ffo 1 + 3 w 


l + 3m 

r 


bj 2 


3w + 5 


2(1 + 3 w) ’ ao-ffo 1 + 


2 l + 3w 

T 3(1+w) 


(58) 


where A , B are constants of integrations and Ji, J 2 denote Bessel functions of the first and second kind, respectively. 
Since we are only interested in the growing mode, we may set B = 0 eliminating the decaying mode contribution. 
Specializing to the case of radiation, w = 1/3, remembering that horizon crossing for the mode under consideration 
(for which k/(aoHo ) = 27t(Ah/Agw) 0 ~ 27t x 10~ 20 holds) happens at thc = 10 40 and using the estimate itsTl) for 
the shear anisotropy at horizon crossing, one determines the remaining constant to be |A| ~ 10 -16 n. Hence, at the 
end of inflation (r = 1), one finally obtains for the sought shear anisotropy 


s » = ©0 ~ 10-45 ' < 59 > 

where the approximation Ji(v,x) ~ x v for small arguments x <C 1 of the Bessel function of the first kind has been 
employed. This is remarkably close to the result one would obtain by simply using the growing mode solution of 
in the limit k/(aoHo) <C 1, that is £(*,) = So T in the case of radiation, which gives So ~ 1CV 46 . If the above value for 
the shear anisotropy is used in car the gravito-magnetic amplification is completely negligible: 


1 

10 



(60) 


We stress that the efficiency of the mechanism depends crucially on the ratio between the coherence length 
Ag of the initial magnetic field and the initial size of the horizon Ah- This ratio, however, disappears when 
the infinite-wavelength limit is taken (see section mvuiD ). Even though the solutions Pi . Pi show a growth 
proportional (quadratic) to the scale factor, the factor of proportionality (cr/H )o (~ 10 -26 or ~ 10 -45 in our first 
and second examples, respectively) is far too small in order to achieve an effective amplification. It follows that 
the interaction between GWs and on average homogeneous magnetic fields is completely negligible in the limit of 
infinitely long-wavelength gravity waves. 
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VI. COMPARISON AND DISCUSSION 


The interaction between GWs and magnetic fields in the cosmological setting has recently been investigated in m. 
where the weak field approximation {Til ] was used. Here one allows for a weak magnetic test field B a in the background, 
whose energy density, anisotropic stress and spatial dependence have negligible impact on the background dynamics: 
B 2 <C /i and Tr a b = —B <a Bb> ~ 0 ~ D a Bb to zero order. In order to isolate linear tensor perturbations, it is necessary 
to impose D a B 2 = 0 = e a fec B b curl B c in addition to the standard constraints u> a = 0 = D a /x = D a p associated with 
pure perfect fluid cosmologies. In the weak field approximation, the main equations governing the induced magnetic 
field arising from the interaction between a weak background magnetic field B a and GWs were derived in 0 for the 
case of a spatially flat Universe with vanishing cosmological constant A and a barotropic equation of state p = wp: 


B , 


W 


-QB, 


W 


1 £ 2 1 

(1 - w) 0 2 + _ 

6 a 


B(i) = 2 


O-(fe) 


-0CT, 


(fc) 


B, 


(«1 


(f) 2 - 


(61) 


where the GWs are determined by the shear wave equation 


a (k) 


-0(71 


(fe) 


|(i-3»)e= + g n 


a (k) 


= 0 


(62) 


Here, the shear is harmonically decomposed as cr a b = &{k)Q a b ’ w hil e f° r the induced magnetic field B a = B^Va 
with Va ] = Q^bQ\ n ) was adopted. The background magnetic field evolves as B a = B°(ao/a) 2 and = ZJjUj Qa' 1 is 
assumed. 

We want to compare our results with the corresponding ones in the weak field approximation. For simplicity, we 
restrict ourselves here to the case of dust. As pointed out above, the only allowed magnetic wavenumber for the 
interacting magnetic field is n = 0, when D 0 Hf, = 0, which leads to £ = k. The published solution for the generated 
magnetic field in the weak field approximation, e.g., equation (21) in 0, however, is not applicable in the limit 
n —> 0. This can be traced back to the choice for the initial conditions for the generated magnetic field made by the 
authors of 0 when solving equations (ED E21, see equation (19) in 0 

In what follows below, we solve equations ED El again, including the full solution for the shear instead of merely 
keeping the dominant part as done in 0. We specify the initial conditions by choosing for every mode k of the 
shear a^)(a o) = oo, <7(fc)(ao) = 0 and for every mode £ = k of the generated magnetic field B^(ao) = 0 = 

Note that this choice of initial conditions differs from that in 0 but agrees with our choice made in section IV. The 
solution, including the background field, for an arbitrary wavenumber k of the shear has the structure 




W 

Dust V 


(a) = B 0 1 + f {Va; k) + O (a ^ 


(63) 


where the function f(^/a;k) is built of several oscillatory terms with amplitude (Agw/Ah)o a t most and the un¬ 
displayed part falls of at least as fast as a -1 / 2 . If this is compared with our result ED , one observes that it differs by 
having another time behaviour. More strikingly, however, is that now the term / k) not only amplifies the seed 
field but also grows like \fa in the long wavelength limit ( k/aoHo <C 1). This is in clear contrast to the gauge-invariant 
result ED , where the seed undergoes amplification but then still decays adiabatically like a 2 . On the other hand, 
in the infinite-wavelength limit ( k —> 0), the exact full solution is now 


B 


(o) 

Dust 


(a) = B 0 



1 + ^ 
H 0 



36 

~5~ 



(64) 


Again, we obtain a solution whose time behaviour differs from that found in ED) - However, the weak field solutions 
agree with our presented solutions in the infinite-wavelength limit when only the dominant part of the solutions is 
considered, at least in the examples considered above. The reason why the solutions obtained within the weak field 
approximation are in general not equivalent to the solutions found using the gauge-invariant approach developed in this 
paper results from the non-gauge-invariance of the weak field approximation, where the magnetic field B a interacting 
with the GW is treated as a weak background field. However, gauge-invariance requires B a to vanish exactly in the 
FLRW background. We remind the reader once more that our procedure solves firstly for the gauge-invariant variable 
/ 3 a = B <a> + |0B a , from which the magnetic field B a measured in the frame of reference of u a can then subsequently 
be found. 

A further important remark concerns the issue of conductivity. We have seen earlier that, within our assumptions 
and for spatially flat Universes, the gravito-magnetic interaction leads to an induced magnetic field which is inde¬ 
pendent of the conductivity of the cosmic medium. This is due to the fact that the interaction does not generate 
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rotational electric field modes which might affect the magnetic field. In the weak field approximation, however, the 
situation is completely different. If one assumes that the conductivity of the cosmic medium that high so that electric 
fields are quickly dissipated away, yielding a curl-free induced magnetic field, then equation m no longer applies 
and one simply has to use 


B {e) + 2 -QB {i) =^ k) B^ (^) 2 (65) 

instead, while the equation for the shear (El is unaltered. This means that the weak field approximation produces 
the same result as our gauge-invariant perturbation approach in the high conductivity limit, and for that case only. It 
is therefore evident that in the weak field approximation the conductivity of the cosmic medium has a crucial bearing 
on the generated magnetic field, in stark contrast to the result of our gauge-invariant approach (see also |28$). 


VII. CONCLUSION 

In this paper we have investigated the properties of magnetic fields in the presence of cosmological gravitational 
waves, using a two parameter perturbation scheme. Using proper second-order gauge-invariant variables (SOGI), we 
were able to obtain results in terms of clearly defined quantities, with no ambiguity concerning the physical validity 
of the variables. The full set of equations determining the evolution of the gravitational waves and the generated 
electromagnetic fields was presented, and the integration shows an amplification of the induced magnetic field due to 
the interaction of a ‘background’ magnetic field with gravitational waves. The magnitude of the original magnetic 
field is amplified by an amount proportional to the magnitude of the gravitational wave induced shear anisotropy 
and the square of the field’s initial co-moving scale Once the amplification saturates, the magnetic field dissipates 
adiabatically as usual. The results were discussed in different fluid regimes, in particular dust and radiation, and 
it was established that the dominant contribution to the magnetic field is the same in both fluid regimes. We find 
that the magnitude of the gravitational boost depends significantly on the manner in which the estimate of the shear 
anisotropy at the end of inflation is calculated. For a seed field of 1CU 34 G spanning a comoving scale of about 10 kpc 
today, the shear anisotropy at the end of inflation (during which we assume H ~ 10 13 GeV) should be larger than 
10 —40 for any noticeable amplification of the seed field to arise at all. 

Moreover, we further recalculated the induced magnetic field employing the weak-field approximation, thereby 
extending previous results in [Tfij . and compared the solutions with ours derived in a gauge-invariant manner using 
SOGI variables. It was found that there is a significant difference in the growth behaviour of the magnetic field when 
SOGI variables are used as compared to the case of a weak-field approximation scheme. While the two methods agree 
in the limit of high conductivity, they seem to be compatible otherwise only in the limit of infinitely long-wavelength 
gravitational waves when the dominant part of the solution is considered. 
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APPENDIX: COMMUTATION RELATIONS 

Here we present various commutator relations which have been used in the text. The relations are given up to second 
order in our perturbation scheme. The vanishing of vorticity, u ab = 0, is assumed throughout in conjunction with the 
constraints D a fx = T) a p = 0 which isolate the pure tensor modes. All appearing tensors are PSTF, S a b = S <a b> , and 
all vectors V a , W„ are purely spatial. 

Commutators for first-order vectors V a : 

(D a V b ) ± = D a V b - i(-)D a Uf, - a a c D c V b + H d e dbc V c (A.l) 

(curl V a ) ± = curl U a - ±0 curl V a - e abc a bd D d V c - H ab V b (A.2) 

D[ a D b ]I4 = [±0 2 - |(^ + A)] V [a h, b]c + (^@(T c[ a- E c[a ) V b] 

+h c [ a ( E b ] d - 3©cTb]d) V d 


(A.3) 
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Commutators for first-order tensors S ab - 

(D a S bc ) ± = D a S bc -±QD a S bc ~a a d D d S bc + 2H a d e de{b S c) e (A.4) 

(D b S ab ) ± = D b S ab — ^QD b S ab — <j bc D c S ab + e abc H b d S cd (A.5) 

(curl S ab ) ± = curl S ab — curl S ab — a e c e cd ( a D e S b ^ d + 3H c<a S b > (A.6) 

curl curl S ab = -D 2 S ab + (n + A - ±0 2 ) S ab + %D <a D c S b>c 

+35' c <q ( E b> c - i0 CTh> c ) (A.7) 

Commutators for second-order vectors W a : 

(D a W b ) ± = D a w b - ±0D a W b (A.8) 

D [a D b] W c = [I0 2 -1( M + A)] W [a h b]c (A.9) 

curl curl W a = -D 2 W a + D 0 (div W) + § (n + A - ±0 2 ) W a (A.10) 
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